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Abstract 

In our previous paper (Eur. Phys. J. E 4, 121 (2001)) we proposed a coarse- 
grained elastic energy for nacre, or stratified structure of hard and soft layers. 
We then analyzed a crack running perpendicular to the layers and suggested 
one possible reason for the enhanced toughness of this substance. In the 
present paper, we consider a crack running parallel to the layers. We propose 
a new term added to the previous elastic energy, which is associated with 
the bending of layers. We show that there are two regimes for the parallel- 
fracture solution of this elastic energy; near the fracture tip the deformation 
field is governed by a parabolic differential equation while the field away from 
the tip follows the usual elliptic equation. Analytical results show that the 
fracture tip is lenticular, as suggested in a paper on a smectic liquid crystal 
(P. G. de Gennes, Europhys. Lett. 13 (8), 709 (1990)). On the contrary, 
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away from the tip, the stress and deformation distribution recover the usual 
singular behaviors (y/x and 1/y/x, respectively, where x is the distance from 
the tip). This indicates there is no enhancement in toughness in the case of 
parallel fracture. 

87.68,+z, 46.50.+a, 83.70.Dk, 81.07.-b 
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I. INTRODUCTION 



A number of living soft matters derive their strength from some composite structures. 
Tooth and timber are among the examples. These composite structures have motivated 
studies mainly in industry or technology-oriented field to produce strong materials, such as 
raw material for construction or air plane and automobile tires. [|T]-0] One of the purposes 
of this paper is to present some complementary understandings of such a problem from a 
viewpoint of a physicist. 

In the first of this series of papers, || which deepened the scaling picture presented in 
f|, we studied nacre which has an alternating laminated structure of hard and soft 
layers (Fig. 1). One of the significant features of this substance lies in the length scale of the 
layers. The thicknesses of hard inorganic layer and that of the soft organic layer d s are of 
the order of nanometer and micron, respectively. This allows us a coarse-grained treatment 
in which average deformation field may be well defined and can be used for the analysis of 
the field even near the tip. Simplifying the coarse-grained elastic energy further in terms of 
Fourier components of the deformation field, which is somewhat different from conventional 
approaches, we could show analytical solutions and thereby demonstrate one possible 
mechanism for the toughness of nacre for a perpendicular crack. 

This paper does not concern the enhancement in toughness as opposed to the title. 
Rather, we proceed to study a case of parallel crack with the fixed grip condition from the 
continuum view (as for non-continuum treatment on laminated composites, see for example, 
Toll ) . As we see later, this coarse-grained treatment again leads to a reasonably simplified 



picture of the problem. To describe the fields near the tip, we shall find that an extra term 
is required in our elastic energy. This term introduces a new length scale A, as we see below, 
and the extra term becomes important near the tip (x <C A, where x is the distance from 
the tip). This term is suggested previously concerning some smectic liquid crystal [11J and 



a lenticular tip form is predicted from scaling arguments. |12j Here, we obtain an analytical 
solution to the stress and deformation distribution near the tip; the tip form in this case 
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takes also the lenticular form. The appearance of a parabolic differential equation for the 
deformation field due to the extra term is another feature of this paper. We also obtain an 
analytical solution for x ^> A, which allows us to estimate the fracture energy; there is no 
energy enhancement within our treatment. 

The Young modulus of the hard and soft layer are denoted E h and E s , respectively, 
where 

E s = eE h . 

In the following, we consider the case of parallel fractures under the plane strain condition, 
i.e. e zx = e zy = e zz = (Fig. 1). Here, is the strain derived from the displacement field 
u i} i.e. 

_1 (dm duA 
13 2\d Xj dxi)' 

where (xi,X2,xs) = (x,y,z). 

The ensuing analysis is based on the conditions appropriate for nacre: 

e < 1, 
d s < dh, 
e = ed/d s 1, 

where, d = d s + d h ~ d h . 



A. Elastic energy for the nacre in the thin layer limit 

Assuming that the thickness of layers are thin in the sense that we can neglect the stress 
change over a few layers, we can introduce macroscopic strain field to have the following 
elastic energy M 

, _ E 2 £-0 2 Ep 2 vEp 

1 ~ 2(1 - u 2 f xx 2 6yy 1 + u 6xy 1 - v txxe ^ [i) 
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where 



E = E h , 
Eq = eE h . 

We have assumed that the Poisson ratio is the same for both layers for simplicity (even if 
£, > or £ < 1, i/ s and are typically of the same order). We have checked (for the 
perpendicular crack) the stress distribution resulting from this energy exhibits small change 
over a few layers, which justifies the above assumption. 

The strain-stress relation results from this energy by the relation o"y = Of ' jdeif 



E v 

@XX Z o&XX H~ Z EoCyy, l^^V 

1 - v 1 1 - V 



V 



Oyy = E [e yy + Y~^ eXX ) ' 

o m = jfL^. (2c) 



B. Parallel fracture: inclusion of the bending effect 

Near the fracture tip, the stress gradient is large; we may need to include the stress 
changes within a single layer in some cases. In the case of parallel fracture, the bending, 
which originates from the non-uniform stress, is important as we see below. The correction 
due to the bending effect is given by (neglecting the bending energy of the soft layer) 



B 



B h ( d 2 u y \ 2 _ K ( d 2 u y \ 2 



2d \dx 2 J 2 V dx 2 
where the bending moduli is given by [|1^ 



Dh 



12(1 -v 2 ) 
The total energy is then given by 

K /<9V N 2 



As we see later, in the case of parallel fracture, the first term in / in Eq. ([I]) becomes 
negligible; the remaining terms in / are all associated with the weak modulus E s while the 
bending term is associated with the strong modulus Eh- This suggests the possibility that 
this higher-order derivative term make a significant contribution to the energy in the parallel 
fracture. 

In accordance with the introduction of the bending term, we have another length scale 
A as announced: 

A 2 = — d l 
E ~ e ' 



II. SCALING PREDICTION 

As we see below, at large scale where A C i, we have u x <C u y , and the elastic energy is 
simplified to 

f= Eo( duy\ 2 E fdu y \ 2 

J " 2 ^ dy J 4(l + i/) \dx J ' { ! 

At small scale A ^> x, we also have u x <C u y and the elastic energy is simplified to 

E (du y \ 2 K (d 2 u y V 

This energy in Eq. fl5|) has been discussed for smectic liquid crystals. [ll],[l2| By consider- 



ing energy balance as in the Griffith theory, |L4| it is predicted that the system exhibits a 



"lenticular fracture," where the opening angle 9q is finite, i.e., 

where Gq(= 27 ) is a thermodynamic separation energy for the soft medium (where 70 is 
the surface energy of the organic substance). Note here this angle is small (8 <C 1) because 
Go/Eq ~ a , where a is a molecular length; thus we see the angle is small {9q ~ a/A). (Here, 
the relation, Gq/E ~ a, results from the following argument. The scaling expression for 



6 



the elastic modulus for polymer melt is given by E ~ cT/N e where c ~ 1/a 3 and N e ~ 100, 
| 15| while the surface tension is typically given by 7 = T/a 2 . [[16] Thus, we have E ~ 7/0.) 
In the inner (small scale) region, the deformation of the free surface is given by 

u s = u y (x,y = 0) = 9 x. 

In the outer (large scale) region where the bending term is neglected, we expect the usual 
singular behavior for the stress and deformation distribution because the deformation field 
satisfies the two-dimensional Laplace equation resulting from the energy in Eq. (f|); 



a s = a yy (x,y = 0) 



u s = —Vx. 

-C/n 



-X 



-0 

Requiring the matching of the deformation field at x = A, we have 

This shows that the stress intensity factor is the same order with that for a pure organic 
system; there is no enhancement of the separation energy. 

Kj usually depends on the remote tensile stress and, at the moment of fracture (where 
Coo = & Failure), the square of this coincides with the elastic constant multiplied by the 
fracture energy (e.g. K] = jiG). In the above, it seems that Kj does not depend on the 
remote stress from the beginning. This is because the expression (|6|) is derived from an 
energy-balance consideration which requires the condition "at the moment of fracture." 

III. SIMPLIFICATION OF THE ENERGY AT LARGE AND SMALL SCALES 

At equilibrium, by requiring 5F = F{u x + 5u x ) — F{u x ) = for any displacement Su x , 
where 

F= fdr( E e 2 +^e 2 + ^e 2 + e + £ X) (7) 

J 2(1 - v*) xx + 2 y y + 1 + v xy + 1 - v xx yy + 2 1 dx 2 J ' K } 
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we have 



where 



a 



2 (1 + u) 

For the displacement of u y , we have 



d 2 u x d 2 u x d 2 u y 

ITT + ea ^TY + ^TTTT = °> 
aar cn/ z oxoy 



1 1 / 1 2z/ 



2 VI + i/ l-i/ 



d 2 u y t n d 2 u y t Q d 2 u x ^ 2 d 4 u y _ 



+ a- 



dy 2 dx 2 dxdy 
Let us try to find a solution of the form 



<9x 4 



u x = uexp(ipx) exp(iqy) 
u y = vexp(ipx) exp(iqy) 



and put them into the above two equations. We have 

/ 



p 2 + eaq 2 



\ 



eppq 

(3pq q 2 + (a + \ 2 p 2 )p 2 
Requiring (u, v) ^ (0, 0), with a' = a 2 — (3 2 , we have 



) 



v / 



eaq 4 + (l + e (a' + a\ 2 p 2 )) p 2 q 2 +(a + A 2 p 2 ) P 4 = 0. 

Since we concern the region, p <C 1/d, for a continuum theory, from the relation A 2 ~ d 2 /e 
we have eA 2 p 2 <C 1 and e\p <C 1. Then, we arrive at the two solutions, 

2 

g 2 = -p 2 ( a + AV) + 0(e), (1 + 0(e)) 

ea 

The minus sign on the right-hand side indicates that the solution is a damping function in 
the y direction and some trigonometric function in the x-direction. 

From the second solution, we have v ~ \feu. Since we seek the solution with v > u in 
our boundary condition, we chose the first, i.e., q 2 = —p 2 (a + \ 2 p 2 ). Thus, at large scale 
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(Xp <C 1), we have q 2 = —p 2 ct, which leads to u ~ et>. These facts allow us to reduce Eq. 
(0) to the form announced in Eq. ([|). At small scale (Xp ^> 1), we have q 2 = —X 2 p A and, 
thus, u ~ eAjw, which allows Eq. (|7|) to be reduced to the form in Eq. (|[). 

Unlike the perpendicular case, [|J the energy is not simplified into a single form over all 
length scale in the present case. Instead, we observe the simplification of the energy at both 
the large and small scales, i.e., Xp <C 1 and Xp ^> 1. We note here that A is rather large 
compared with the layer thickness, i.e., A ~ d/^/e. The ensuing analysis is based on the 
condition, tf< A<L, where 2L is the y dimension of the sample. 



IV. LARGE SCALE SOLUTION 



At equilibrium, by minimizing the energy given in Eq. (|j), we have 

( d 2 d 2 \ 

+ TT^T ]U V = 0, 



2 T o 2 "2/ 



dx 2 dy 



where 



x 



y/2(l + u)x. 



When the crack tip is located at the origin of the x-axis the appropriate boundary conditions 
are, under the fixed condition, as follows: 

u y = ±u at y = ±L, 

u y = for y = 0, x < 0(x < -A), 
Qu 

— - = for y = 0, x > 0(x > A). 

Note here that these boundary conditions specify the conditions at the points away from the 
origin. As in the perpendicular case, in analogy with the boundary problem for a variable 
condenser, 0] we have the displacement; 

2u 



u y = Im 



7T 

with 



log ^ 2 /(2L) + ^ 2 A_l) 1/2 ^ ; (8 ) 



z = y + i (x - x ) . 

In the above, the branch of the function z 1 ! 2 is chosen such that z 1 ^ 2 = r 1 ! 2 exp(i9/2) with 
< 9 < 2tc (e.g., (— l) 1 / 2 = i). Here, the remote boundary conditions allow the shift, 
(= xq/ \J2{1 + u)^j , in our solution, as far as x <C L. 

In the following, we consider only the region y > 0, since the problem is symmetric with 
respect to the x-axis. When \z\ <C L, we have 



u y (x,y) = — Im 

71 



(i?rz\ 

\~TJ 



1/2 



(9) 



In the vicinity of the origin (but still x >> A), it reduces to a parabolic form: 



u y (x,y = 0) = 2m 6 ) (x)i 



'rr — x 
ttL 



(10) 



where is the Heaviside step function. Here and hereafter, the notation y = should be 
understood as the limit where the variable y approaches to an positive infinitesimal quantity. 
This is because the consideration for the region y > is enough due to the symmetry. 
The nonzero component of the stress field is given from a yy = E e yy ; 

e inz/(2L) 



®yy &oc 



(ii) 



where the remote tensile stress is given by 



771 «0 



When \z\ <C L, we have 



a yy (x,y) = doc Re 



and stress distribution near the tip is given by 

o- yy (x,y = 0) = Ooo^-rr) 



t) 



-1/2 



~(X- X )) IT 

Note here that these expressions are valid only when x > A(~ d/y/e) where \ <^ L. 



(12) 
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V. SMALL SCALE SOLUTION 



From Eq. (j^), we have a "pseudo bi- diffusion" equation, 

?y = A 2 ^. (13) 

dy 2 dx A 

The solution to our problem has to satisfy the following two boundary conditions at least 

u y = for y = 0, x < 0(\x\ <C A), 
du 

—1 = o for y = 0, x > 0(x < A). 
In addition, the solution should be able to match with the large scale solution. 

A. General solution for the boundary problem 

The field u y satisfies the (pseudo) bi-diffusion equation, 

(dy + Xdl) (dy - Xdl) Uy = 0. 

Putting S u = (d y — A<9^) u y , we see that S u is a solution to the anti-diffusion equation 

(dy + Xdl) S U = 

and that u y satisfies the inhomogeneous diffusion equation 

(d y - Xd^ u v = S u - 

The source function S u is a solution of the anti-diffusion equation (a typical solution is 
exp (^j) I s/^irXy) and blows up when y approaches to zero, except that S u is independent 
of y, i.e., S u = A + Bx (this fact can be confirmed later) where A and B are constants 
independent of x and y. 

Since u y satisfies the bi-diffusion equation, o~ yy also satisfies the bi-diffusion equation. 
Putting So- = (d y — Xdl) a yy> we have 



(d y - Xdl) °yy = Sa- 



il 



To have a physical solution, we have S u = A + Bx and we conclude that S c = because 
S a = E d y S u . Namely, a yy is the solution of the diffusion equation. As we see later, the 
function S u obtained starting from the diffusion equation for a yy is actually in the form, 
S u = A + Bx. In this way, we shall obtain the general solution for u y . 

The general solution for the initial boundary condition (here, "initial" corresponds to 
y = 0) 

a yy (x,y = 0) = F(x) (14) 

is given as 

\2\ 



(x — x'Y 
y/4nXy J~oo v ' v \ AXy 



a yy (x,y) = jz z J F(x')exp ( — ) dx' . (15) 



In the following, we determine the source function F(x) considering appropriate boundary 
conditions. Requiring the boundary condition 

o~ yy (x, y) = for y = and x > 0, (16) 

we may put 

F( x ) = 9{-x)f(x) + a\5(x) (17) 

where 5(x) is the Dirac's delta function and f(x) and a are the quantities to be determined. 
As we see later, the second singular term in the source function is required to have a smooth 
matching to the outer solution. In this way, we have 



a ^ y) = vmL f{x ' )e ^ ^^i 1 -) dx ' +a Hk exi> (-4^) • (18) 

By integrating Eq. (|I~8|) over y, we have (erfc (x) — 1 — erf(x), erf(x) = Jq e~ tZ dt) 

„ f° r, r\ f I V (x-x'f x — x' ( x — x'\\ . 

E u v ( X , v )=J_J( X )^-e «. -—^[-^jjdx 



+ - [ ^e,p (-^-) - ferfc ) | + (19) 
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where we have used the formula given in |17|,18|, in which formulae the integration constants 



are chosen so that the result vanishes at x = — oo. In the above, C(x) is the solution of the 
bi-diffusion equation independent of y, i.e., 

C(x) = C + C\X + c 2 x 2 + c 3 x 3 . 

The boundary condition required for u y (x,y) is 

u y (x, y) = for y = and x < 0. (20) 

To determine the function f(x) to accommodate this boundary condition, we consider 
u y (x,y) at y = 0. Noting the relation, 



dx + / dx' J f(x') — erfc 

-00 Jx / ZA 



where the first term in the right-hand side is zero at y = and the second term is nonzero 
only when x < 0, we have 

E u y (x, y = 0) = 9 (-x) Q £ f(x') (x - x') dx' - ax^j + C(x). 

Here, note the relation erfc ( (x — x') / y/AXyj — ► 26 (—x) when y — > + . 

Since 6{—x) = 1 — 9(x), the boundary condition (p0|) is satisfied only when 

C(x) =ax- \ T f(x') {x - x') dx', (21) 



A Jo 

to have 



E Q u y (x, y = 0) = 9 (x) (ax — — J f{x) (x — x) dx^j . 

The condition is strong in the sense that C(x) is the polynomial of x with the order 
lower than the fourth and is independent of y; f(x) is at most linear in x, 

f(x) = do + a\x. (22) 
13 



In this way, we obtain 

E u y (x, y = 0) = 9 (x) (ax - - Qa x 2 + jUi^ 3 )) ( 23 ) 
and from Eqs. ( |I4D and ([T7|), we have (for x ^ 0) 

ayy(x,y = 0) = 0(— z)(a + aisc) (24) 

VI. MATCHING BETWEEN INNER AND OUTER SOLUTIONS 



Eq. (|23| ) contains the three unknown constants, Oq, Oj, and a, while the outer solution the 
one unknown, xq~, there are four unknowns. To determine these constants we require that 
u y and d x u y match the outer solution at x = A while a yy and d x a yy match the large-scale 
solution at x = —A. 

A. Outer solution 

Introducing a (which is different from a used in Sec. Ill) 



a 2 = J2(l + u) 



we have 



u y (x,y = 0) = 2u a9(x)] 



IX — Xq 

ttL 



S-Uv(x,y = 0) = 9(x) u ° a 

a/7tL [X — Xq) 



U0 

^0 ' a^ixL (- (x - Xq)) 



-<jyy(x,y = 0) = 9(-x 
^r^-a yy (x, y = 0) = 9(-x) U ° (- (x - x )y 3/2 
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B. Inner solution 



Introducing a, b and u (which is different from u in Sec. Ill) 

Aa , a\\ 2 a\ 
a = — — , o = — — , u = — 

-EqM -KqW Eq 



we have 



y = 0) = u6{x) 



x I a 2 ^ 3 



d ( 1 ( a b 2 

^(^ = 0) =«*(*) + — x 

■^r<r m (x, y = 0) = + ^ 



C. Matching 

We have four conditions, while we have four unknowns, a,b,x , and u, for the four 
conditions. 



At x = A 



Here, we have introduced 



At x = -A : 



fa b \ 1 aMol 
" M U + 2aJ +M A = V^ 



/ / A i x ° 



It , , W 

(a — o) 



A aAvY 
6 



2aAA 2 (v^ 7 )' 
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Here, we have introduced 



7 =1 + T , 



from which we have — x + x = A(l + x /X) = Xj'. With the definition, 

77 = — I, 
u 

we have 

67-I 



a = 4 — 2a?7- 

4^-1 

6 = -6 + 6ary - (25) 



and 



2 7 ' + 1 
a = 77 3- 

2a(V)* 
6 " 



2a(V) 3 

Thus, we have a set of equations for 7 and 77 : 

67-I 2V + 1 
4 — 2ar] — — = r] — _ n a 



V7 2a(V) 5 
7- 1 = g 
\/7 2a(V)i 



-6 + 6ar/ 47 1 = — ^—3- (26) 



Numerical solution to this set of equation with \x \ < X (or |1 — 7I < 1) is 7 = 
0.71842, 77 = 0.38752 for 1/ = and r] = 0.378 70, 7 = 0.6656 for v = 1/2. 

Approximate but analytical expressions for them can be obtained as follows. Introducing 
77' = I/77 and 5 by 

5 = x /A,7 = 1 - 5, V = 1 + 5 

and, assuming, 5< 1, the set of equation can be reduced to 

7 A 6-55 



477' - 2a (5 - -5) = 
-677' + 6a ^3 - ^5 



4a 

5 \ 2-35 



4a 
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Solving this equation we have an approximate solution: 

2 (12a 2 - 11) 
~ 3 (12a 2 - 7) 
, _ 48a 4 - 8a 2 - 1 
77 ~ 3 (12a 2 - 7) a 

For a = ^/2(1 + v) with < v < 1/2, 5 ~ 0.4 which barely justifies the assumption, i5< 1, 
while rf ~ 2.4. In this way, we confirm that 7 and t] are both positive quantity of the order 
of unity for any value of v (0 < v < 1/2). 



VII. OVERALL SOLUTION 



For \x\ < A, from Eqs. fll9D , fl21|) , and ( p2|) we have 

^-Ki- 2+ H( i+erf (v§)) (27) 

+ l -Y (a + bX) erfc ( *) — (-bY + + ,/Z e -£ 



2 VV4F/ \3 X 7 V 7T 

and from Eqs. (|i~8|) and (E3) we have 



o- yy (x,y)\ I. ( X \ Y _xi 

— = - (a + 6X) erfc — == - b\ — e «- 

where a (~ 0.4) and b (~ 0.1) are given by Eq. ( P5| ) and w = MoV 7 ? (~ 2.6m /) with 7 and 77 
given by Eq. fl26|) . Here, X = x/X and Y = y/X. 

For |x| > A, w y and cr w are given by Eqs. © and ([□]) where a; = (1 — 7)A (~ 0.3A). 

As expected and clear from (|27|), these expression are only valid for y -C A (at small 
scale, we have q 2 = —X 2 p 4 with Ap ^> 1, from which we have q 2 3> p 2 ^> 1/A). 

If we operate (<9 y — AS 2 ) to the right-hand side of Eq. ( p7|) we obtain a + bX, which 
corresponds to S u in the previous section and thus confirm the statement mentioned there. 

A typical shape of the crack tip and stress distribution around the surface (0 < y <C A) are 
given in the plots in Fig. 2 with the parameters given in Discussion. The solution becomes 
precise only when i C A or a; > A. Near the tip (x <C A), the tip is lenticular (~ x) 
while, away from the tip (x ^> A), the deformation field recovers the ordinary square-root 



law (~ y/x) . ||-gT 
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VIII. DISCUSSION 



Since we have the usual singular behavior for the stress and stress fields we have 

G ~ (™)large scale ~ a ooU . (28) 

On the other hand, from the Griffith's energy balance, we have 70 ~ Lc^/Eq. [[J Thus, the 
fracture energy is the order of 70 as predicted by the scaling consideration. 
Since we have 

M x ) = -§^ u y( x >y = 0) = u9(x) ^ - (j^x+ 2X3^)) 
the crack tip angle at x = is 

e„ s 5 o( o) = u -l ~ 

By rewriting the last expression with the aide of (p8|), we have 



' 7o 
E X 



0o ' 

as predicted in Eq. (|). 

Typical parameters for nacre are as follows: 

d = 1/zm 
L = 1 cm 
e = 1/250 

For this parameter we have 

A = 250/im 
L/X = 40 
d/X = 1/250 

These numerical values allow us the continuum description with length scales A and L since 
the relation, d <C A <C L, holds for these values. 



If we did not include the last term in Eq. ( |I7|) , we would not have a linear term in the 
deformation field and it would contradict with the scaling expectation (in addition, we could 
not make a smooth matching as we see below); we would have 

1 u ( b \ „. , , 

—a yy = j\a + -x\ 6{-x) (29) 

Here, we have only two unknowns; a 1 = ua and b' = ub. Let us assume first that a' ^ 0. 
Then, the first equation implies a' < for a positive u y around the origin (x = 0), while 
the second equation implies a' > for a positive a yy (We have assumed u > 0). Thus, for 
consistency, we have a' = 0. 
Thus, we have 

u y (x,y) = -—x 3 9(x) 

For positive u y and <r w , we have b' < 0. By introducing the shift Xo for the outer solution 
as before, the matching of the both fields themselves is possible. However, since the inner 
solution u y (~ x 3 ) is concave (positive curvature) while the outer solution (~ y/x) convex 
(negative curvature), the smooth matching is impossible. For a similar reason, the smooth 
matching of the stress is also impossible and the stress field at x = tends to zero. These 
facts suggest the need for the last term in Eq. ([17|) . 



IX. CONCLUSION 

In this paper, we have developed a continuum theory for a nanoscale layered structure 
with a parallel crack. We indicated the importance of the bending term near the crack, which 
leads to a parabolic differential equation and, thereby, a lenticular tip form. From the coarse- 
grained view, we have a simplified picture with analytical expressions; the deformation field 
is linear function of x (the distance from the tip) near the tip and goes back to the usual 
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square-root function of x away from the tip. As a result, the fracture energy of nacre against 
the parallel crack is of the same order as the soft organic material; there is no enhancement 
in toughness within our simplified view. 
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FIGURES 

FIG. 1. Nacre structure: the (inorganic) hard layer thickness dh is of order of micrometer while 
the (organic) soft layer thickness d s nanometer. The y-axis is perpendicular to layers and the 
sample is long in the z-direction. The cracks in the y — z plane and in the x — z plane are called 
the perpendicular and the parallel fractures, respectively. 

FIG. 2. Crack shape (a) and stress distribution (b) for a parallel crack. The bold line corre- 
sponds to the overall solution while the fine line and the broken line correspond to the large-scale 
solution and the small-scale solution, respectively. The overall solution becomes exact only when 
x < A or i > A. Near the tip (x <C A), the tip is a linear function (~ x) while, away from the tip 
(x » A), the deformation field recovers the ordinary square-root law (~ x 1 / 2 ) . The parameters 
used for the numerical calculation are d = 1/xm, L = 1cm, and A = 250/um where e = 1/250, 
L/X = 40, and d/X = 1/250 (See Discussion) 
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This figure "Fig 1. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0203239vl 



This figure "Fig2.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0203239vl 



